Army  Research  Laboratory 


The  Effect  of  Helicopter 
Main  Rotor  Blade  Damage  on  the 
Rotor  Disk  (Whole  Rotor)  Motion 


by  Joseph  Fries 


ARL-TR-2241 


June  2000 


20000614  019 


Approved  for  public  release;  distribution  is  unlimited. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other 
authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  die  use  thereof. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return 
it  to  die  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5068 

ARL-TR-2241  June  20°0 

The  Effect  of  Helicopter 
Main  Rotor  Blade  Damage  on  the 
Rotor  Disk  (Whole  Rotor)  Motion 

Joseph  Fries 

Survivability/Lethality  Analysis  Directorate,  ARL 


Approved  for  public  release;  distribution  is  unlimited. 


Abstract 


When  a  helicopter  main  rotor  blade  is  ballistically  damaged,  an  imbalance  is  created  in  the 
rotor,  causing  the  rotor  disk  to  execute  unwanted  motions,  which  are  detrimental  to  performance. 
The  normally  smooth-flying  helicopter  develops  new  vibrations  that  can  be  physiologically 
annoying  or  debilitating  to  the  pilot,  can  exceed  structural  fatigue  endurance  limits,  can  cause 
aeromechanical  instabilities,  and  can  reduce  helicopter  performance  ability. 

This  report  examines  the  effect  of  the  loss  of  the  outboard  section  of  one  rotating  blade  of 
a  rotor  set  of  four  blades  on  the  fixed-system  (nonrotating)  rotor  disk  motion.  The  report  shows, 
beginning  with  the  rotor  blade  forcing,  how  a  damaged  blade’s  response  changes,  and  how  this 
change  feeds  into  the  rotor’s  fixed-system  disk  motion  (the  disk  referring  to  the  blades  acting  in 
concert  as  a  whole  entity). 

With  a  normally  undamaged  rotor  (referring  to  all  the  blades),  there  exists  within  the  rotor 
itself  the  capability  of  motion  canceling  of  certain  frequencies  depending  on  the  number  of  rotor 
blades  in  the  rotor.  This  study  tracks  each  individual  harmonic  (integer  multiples  of  the  rotor 
speed)  frequency,  one  at  a  time,  in  order  to  obtain  a  first-principles  understanding  of  the 
phenomena  involved  with  rotor  imbalance. 
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1.  Introduction 


When  a  helicopter  main  rotor  blade  is  ballistically  damaged,  an  imbalance  is  created  in  the  rotor, 
causing  the  rotor  disk  to  execute  unwanted  motions,  which  are  detrimental  to  performance.  The 
normally  smooth-flying  helicopter  develops  vibrations  that  can  be  physiologically  annoying  or 
debilitating  to  the  pilot,  can  exceed  structural  fatigue  endurance  limits,  can  cause  aeromechanical 
instabilities,  and  can  reduce  helicopter  performance  ability. 

This  report  examines  the  effect  of  the  loss  of  the  outboard  section  of  one  rotating  blade  of  a  rotor 
set  of  four  blades  on  the  fixed-system  (nonrotating)  rotor  disk  motion.  The  report  shows,  beginning 
with  the  rotor  blade  forcing,  how  a  damaged  blade’s  response  changes,  and  how  this  change  feeds 
into  the  rotor’s  fixed-system  disk  motion  (the  disk  referring  to  the  blades  acting  in  concert  as  a  whole 
entity). 

With  a  normally  undamaged  rotor  (referring  to  all  the  blades),  there  exists  within  the  rotor  the 
capability  of  motion  and  force  canceling  of  certain  frequencies  depending  on  the  number  of  blades 
in  the  rotor.  This  study  tracks  each  individual  harmonic  (integer  multiples  of  the  rotor  speed) 
frequency,  one  at  a  time,  in  order  to  obtain  a  first-principles  understanding  of  the  phenomena 
involved  with  rotor  imbalance. 

It  is  the  purpose  of  this  work  to  “break”  the  damaged  rotor  blade  imbalance  phenomenon  into 
its  most  fundamental  elements,  these  elements  being  the  relative  individual  blade  responses  between 
the  damaged  and  undamaged  blades,  the  blade  (rotating)  forcing  frequencies,  the  fixed-system 
(nonrotating)  rotor  disk  frequencies  and  phase  angles,  the  type  and  location  of  rotor  blade  damage, 
the  number  of  blades  that  compose  the  rotor,  and  the  particular  blade  of  the  rotor  set  that  is  damaged. 
Among  these  variables  are  deterministic  coupled  relationships  that  work  together  to  produce 
fixed-system  rotor  disk  motions  as  a  function  of  the  individual  rotor  blade  response  characteristics. 


Note:  A  list  of  terms  used  in  the  equations  that  follow  is  included  at  the  end  of  this  report 
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Each  rotor  blade  is  excited  by  a  common  periodic  azimuthal  harmonic  forcing  function  (i.e.,  each 
blade  of  the  rotor  experiences  the  same  forcing  at  the  same  position  in  the  azimuth).  However,  for 
a  damaged  blade,  because  of  loss  of  mass,  blade  length,  and/or  stiffness,  its  dynamic  response  is 
different  from  those  of  the  undamaged  blades.  This  difference  nullifies  the  vibration  cancellation 
effect  that  exists  when  all  the  blades  are  similar  or  undamaged.  As  a  result,  these  uncanceled 
motions  at  particular  frequencies  move  into  the  fixed  system  as  new,  unwanted  forcing  of  the 
fixed-system  rotor  disk  motion  that  ultimately  affects  the  aircraft. 

The  present  analysis  uses  a  Newtonian  approach  (Meyer  1969)  in  deriving  the  blade  equations 
of  motion,  with  2  degrees  of  freedom  for  each  blade  of  the  rotor  set.  A  representative  four-blade 
rotor  is  modeled  in  the  analysis.  Each  blade  is  represented  as  a  flapping  mass  on  a  blade  length  with 
a  root  effective  flap  hinge  and  a  root  angular  spring.  In  addition,  a  hub  mass  is  represented  at  the 
flap  hinge  on  top  of  a  vertical  spring.  Each  of  the  blades  is  forced  with  harmonic  vertical  shears  at 
the  blade  tip. 

Since  the  object  of  this  analysis  is  to  study  how  and  which  frequencies  are  transferred  to  the 
rotor’s  fixed  system  when  a  blade  is  damaged,  the  magnitudes  of  the  forcing  function  can  be  made 
arbitrary,  and  only  the  relative  change  in  magnitude  between  the  damaged  and  undamaged  case  is 
important  Therefore,  this  analysis  normalizes  both  the  blade  responses  and  the  rotor  disk  responses 
by  the  undamaged  rotor  blade  response  for  each  individual  harmonic.  However,  it  is  recognized  that 
there  is  a  natural  amplitude  relationship  among  the  harmonics  for  a  blade’s  aerodynamic  forcing. 
The  present  analysis  concentrates  on  understanding  rotor  imbalances  and  vibration  propagation 
fundamentals,  rather  than  attempting  to  quantify  rotor  motions  and  loads. 

In  this  analysis,  trigonometric  products  are  involved  in  the  mathematics  of  transformation  from 
rotating  motion  into  nonrotating  motion  that  results  in  frequency  changes.  For  instance, 
trigonometric  terms,  such  as 


sin  k\|r  cosily 
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occur  when  a  harmonic  such  as  sin  k  \|r  in  the  rotating  system  is  transformed  into  the  fixed  system. 
Here, 


are  different  harmonics,  and 


k  and  1 


\|r  =  £2t 


Q 

is  the  rotor  blade’s  azimuth  position.  Here  \|f  =  Qt,  and  when  Q  is  normalized  by  itself  —  =  1 ,  we 

call  this  1/rev.  For  the  kth  harmonic,  k£2  normalized  again  by  Q  yields  — — ,  which  we  call  k/rev. 
By  means  of  trigonometric  identities  (Riddle  1974),  we  know  that 


and 


sin(k\|r±l\}f)  =  sink\|rcosl\|/±cosk\|rsinl\|/ 


sink\jrcosl\|/‘  =  — (sin  (k  + 1)  \|/ ±  sin  (k-l)Y) 
2 


A  frequency  change  occurs  going  from  the  rotating  blade  system  into  the  rotor  disk  fixed  system. 
We  see  the  rotating  sin  k  \|r  is  transformed  into  two  different  fixed  system  frequency  terms  sin 
(k  +  1)  \|f  and  sin  (k  -  1)  Vf.  This  is  one  phenomenon;  another  is  the  phenomenon  of  cancellation 
when  the  rotor  blades  are  summed  to  give  a  total  rotor  nullifying  effect.  For  instance,  we  may  have 
a  fixed-system  mathematical  form  like 

X  3;  sin  k 

i  =  l 

where 


\\r  +  2n 


Pi  =  flap  angle  of  individual  blade,  and 
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k  =  fixed-system  harmonic. 


If  all  the  blades  are  undamaged  and  physically  equivalent,  the  flap  angle  responses  will  be  identical. 
In  this  case,  it  is  known  that  for  a  four-bladed  rotor,  the  summation  will  sum  to  zero  (Johnson  1980) 
for  harmonics  1, 2,  and  3.  This  is  the  cancellation  effect.  However,  if  there  is  blade  damage  and  all 
the  flap  angle  responses  are  not  equal,  then  no  cancellation  occurs,  and  the  1, 2,  and  3/rev  frequency 
vibration  effects  are  felt  in  the  fixed  system.  These  types  of  considerations  are  examined  in  this 
work. 

The  physical  properties  of  a  generic  (16,000-lb  gross  weight)  single  main  rotor  helicopter  are 
represented  in  this  analysis.  The  blade  weight,  effective  flap  hinge  offset,  hub  weight,  and  the 
associated  springs  yield  an  uncoupled  first  natural  blade  flap  frequency  of  1.05/rev  for  each  of  the 
four  undamaged  rotor  blades. 

Damage  is  assumed  to  be  imposed  on  blade  no.  1  by  removing  the  outboard  portion  (up  to  30%) 
of  the  blade,  thus  resulting  in  blade  mass  and  length  loss.  Both  of  these  losses  cause  the  natural 
frequency  of  the  damaged  blade  to  increase.  This  reduces  the  blade’s  response,  especially  to  the 
1/rev  forcing  from  a  resonance  perspective,  which  causes  the  rotor  imbalances.  There  are  also 
differences  of  response  to  the  other  forcing  frequencies.  These  differences  nullify  the  cancellation 
effect  of  the  normally  undamaged  rotor. 

The  ramification  of  blade  damage  to  the  fixed-system  rotor  disk  tilt  is  that  the  disk  thrust  vector 
is  tilted  in  undesirable  directions  with  undesirable  harmonics.  The  helicopter  trim  condition  when 
the  blades  are  undamaged  is  determined  by  the  steady  tilt  of  the  rotor  disk  in  the  fixed  system  from 
1/rev  blade  flapping  in  the  rotating  system.  A  damaged  blade’s  1/rev  response  is  different  from  the 
undamaged  blades,  and  thus  affects  the  trim  setting.  Beside  the  rotating  1/rev  frequency,  other 
frequencies,  due  to  the  damaged  blade,  get  into  the  fixed-system  rotor  disk  motions  as  harmonics, 
tilting  the  disk  and  the  thrust  vector  back  and  forth  at  these  frequencies,  and  shaking  the  helicopter 
as  vibrations.  There  will  be  physiological  effects  on  humans  subjected  to  these  vibrations,  natural 
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frequencies  of  the  helicopter  will  change  with  possible  instabilities  occurring,  structural  fatigue  life 
may  be  exceeded  by  increased  vibrations,  and  general  helicopter  performance  will  be  degraded. 

Because  of  these  effects,  an  understanding  of  the  underlying  fundamentals  of  rotor  imbalance 
is  important  in  order  to  know  how  to  analyze  their  effects. 

2.  Equations  of  Motion 

The  model  used  for  an  individual  rotor  blade  is  illustrated  in  Figure  1 . 


Figure  1.  Model  of  Single  Rotor  Blade. 


From  Figure  1,  let 

Q  =  the  rotor  speed  (om), 
z  =  vertical  deflection  of  the  hub  mass, 

P  =  blade  natural  angular  flap  displacement  (B), 
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Mi  =  rotor  blade  mass  (Ml), 
M2  =  hub  mass  (M2), 

Kz  =  hub  spring  (Kz), 

Kp  =  blade  spring  (Kb), 
e  =  blade  flap  hinge  offset, 
r  =  blade  radius,  and 
F  =  blade  forcing. 


From  Figure  1,  om  is  the  rotor  speed,  denoted  by  an  end  of  an  arrow  going  into  the  paper  and  a 
dot  coming  out  of  the  paper.  Each  of  the  four  rotor  blades  is  represented  as  per  Figure  1 .  The  final 
equations  of  motion  are 


M^2  -  Mtr 

+ 

"Mt  Q2  r  ((e  +  r)  +  Kp  r2 )  0  ' 

"P' 

•Tl 

►i 

_ i 

-Mjr  Mi  +M2 

y_ 

0  k2 

z 

L  F  J 

where  the  primes  are  the  second  time  derivative  of  |3  and  z. 
Let  the  forcing  function  Fi  of  the  ith  blade  be  of  the  form 


Fj  =  A;  sin  k 


(  Tt  ^ 

\ r  +  2— (i  — 1) 
N 


where 


Ai  =  amplitude, 
k  =  harmonic  number, 

\|t  =  azimuth  of  blade  1 ,  and 
N  =  number  of  blades  in  rotor. 


6 


Solving  for  the  steady-state  solution  of  the  equations  of  motion,  we  have  for  the  first  element  of  the 
solution  of  the  form 

Pi=Ci  sinker  + 

Now,  sum  the  blades  together  and  transform  them  into  the  fixed  system. 

Fixed-System  Lateral  Disk  Plane  Tilt: 

4  (  fi-11  ")  (  (i-1) 

R  =  X  C  sink  \}r  +  2n- -  sin  \\r  +  2n— — 

i=i  l  4  )  \  4 

and 

Fixed-System  Longitudinal  Disk  Plane  Tilt 

4  (  Ci-l)  ^  (  (i-1) 

R  =  I  Cj  sink  v  +  27t^—^  cos  \| r  +  2n±—^ 

y  i=i  ^  4  )  {  4 

Call  the  lateral  fixed-system  disk  plane  tilt  px  and  the  longitudinal  py. 

Since  \|/  =  Qt,  where  t  =  time,  we  need  to  determine  the  Fourier  coefficients  of  both  the  lateral 
and  longitudinal  fixed-system  tilts  to  determine  the  amplitudes  and  frequencies  involved. 

3.  Computer  Software  Analysis  Implementation 

The  implementation  of  the  solution  of  this  study  is  done  in  integrated  steps.  The  rotor  blade 
dynamic  response  solution  is  made  with  MACSYMA.*  This  computer-developed  code  is  contained 

*  MACSYMA  is  an  interactive  system  for  doing  mathematical  computation.  It  handles  numeric,  graphic,  and  symbolic 
calmlatinns  and  incorporates  high-level  programming  language,  which  allows  the  user  to  define  his  own  procedures. 
MACSYMA  Inc.,  20  Academy  Street,  Arlington,  MA  02174-6436.  The  telephone  number  is  1-800-MACSYMA 
(1-800-622-7962,  toll  free  in  U.S.  only);  1-617-646-4550;  FAX  1-617-646-3161;  and  the  email  address  is 
info@macsyma.com. 
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in  a  file  named  dynl.sav.  The  MACSYMA  function  odematsys  is  used  to  solve  for  the  steady-state 
blade  response  solution. 

These  responses  are  used  as  inputs  into  another  program  in  FORTRAN  called  sum.f,  which  sums 
the  individual  blade  responses,  makes  a  coordinate  transformation1  into  the  fixed  system,  and 
performs  a  Fourier  analysis  to  determine  the  amplitudes,  Cj,  and  frequencies,  k. 

4.  Results 


The  individual  blade  harmonic  responses  were  solved  using  MACSYMA,  and  the  results  given 
in  Table  1  were  normalized  by  the  magnitudes  of  the  undamaged  blade  responses  and  are  presented 
in  the  following  tabular  form.  These  are  the  rotor  blade  flap  responses  due  to  an  applied  sine  forcing 
functions  at  the  blade  tips  of  an  arbitrary  amplitude.  The  minus  value  of  - 1.0  in  column  2  for  the 
first  harmonic  is  due  to  the  1/rev  forcing  to  be  below  the  blade  natural  frequency  of  1.05/rev  (i.e.,  it 
is  a  dynamic  resonance  phenomenon). 


Table  1.  Normalized  Blade  Flap  Data  (Rotating  System) 


Sin  Harmonic 

1 

-1.0 

-0.7929 

-0.6151 

-0.4636 

2 

1.0 

0.9064 

0.8129 

0.7182 

3 

1.0 

0.9024 

0.8079 

0.7103 

4 

1.0 

0.9017 

0.8034 

0.7052 

5 

1.0 

0.8971 

0.8074 

0.7112 

10%  loss  of  blade  mass  and  length  at  tip. 


From  Table  1,  the  first  column  is  the  harmonic  of  the  sin  forcing  function  that  is  applied  to  each 
blade  in  the  rotating  system  of  the  set  of  four  blades  from  1  to  5/rev.  Column  2  is  the  blade  response 
for  the  undamaged  condition  normalized  by  the  magnitude  of  the  undamaged  blade  response.  The 
remaining  columns  are  the  normalized  damaged  blade  responses  for  10%,  20%,  and  30%  damage. 
For  instance,  30%  means  a  loss  of  30%  blade  mass  and  30%  of  the  blade  tilt  radius. 
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Next,  the  total  rotor  disk  fixed-system  motion  is  calculated  by  summing  the  contribution  of  each 
individual  blade  and  projecting  it  onto  the  fixed-system  x  and  y  axes.  The  angular  motion  projection 
about  the  x  axis  is  called  BX  and  about  the  y  axis  BY.  Figure  2  illustrates  this  convention. 


Figure  2.  Fixed-System  Rotor  Disk  Tilt  Coordinate  System. 

Figure  2  shows  the  direction  of  the  disk  plane  lateral  tilt  BX  and  the  longitudinal  tilt  BY.  Since 
BX  and  BY  are  time  histories,  a  Fourier  series  analysis  is  performed  on  them  to  obtain  the  sine  and 
cosine  Fourier  coefficients  bk  and  ak,  where  k  denotes  the  harmonic  number. 

The  computed  fixed-system  rotor  disk  tilt  information  is  given  in  Table  2  (undamaged  rotor)  and 
Table  3  (damaged  rotor).  The  results  in  these  tables  indicate  some  interesting  phenomena.  Table  3 
contains  more  entries  than  Table  2;  this  indicates  that  more  frequencies  are  transmitted  into  the  fixed 
system  when  the  blades  are  damaged  than  when  they  are  undamaged.  Thus,  in  turn,  the  rotor  disk 
plane  is  tilting  at  these  additional  frequencies.  The  steady  rotor  disk  tilt  of  -2.0  in  Table  2  has 
changed  to  -1.731  in  Table  3;  this  indicates  that  the  damaged  blade  has  caused  the  helicopter  to  go 
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Table  2.  Fixed-System  Rotor  Disk  Tilt,  No  Blade  Damage 


Table  3.  Fixed-System  Rotor  Disk  Tilt,  Blade  Damage  (30%  of  Tip  Removed) 
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out  of  its  trimmed  condition.  For  instance,  for  the  helicopter  to  fly  in  a  straight,  level  flight  condition 
with  undamaged  blades,  assume  the  lateral  disk  plane  tilt  is  normally  -2.0.  This  value  (-2.0)  is  a 
steady  lateral  fixed-system  disk  tilt  for  die  trimmed  condition.  When  the  blade  is  damaged,  the  tilt 
changes  to  - 1.731,  which  causes  the  helicopter  not  to  fly  in  a  straight  line. 

The  ambient  4/rev  rotor  disk  tilt,  as  seen  in  Table  2,  for  the  undamaged  blade  with  a  value  of  2.0 
or  -2.0  is  a  tilting  that  is  normally  present  in  the  rotor  disk.  With  the  damaged  blade,  the  4/rev  tilt 
is  still  present  but  is  reduced  to  1.85,  as  seen  in  Table  3.  This  means  that  the  4/rev  vibrations 
transmitted  into  the  fuselage  are  reduced,  which  is  a  benefit,  but  at  the  detrimental  cost  of  other 
frequencies  (1, 2,  3, 5, 6,  etc.)  appearing. 

In  Table  2,  the  numbers  across  the  top  are  the  harmonics  in  the  rotating  system,  while  the 
numbers  in  the  left-hand  column  are  the  harmonics  transmitted  into  the  fixed  system.  For 
undamaged  blades,  we  see  that  only  the  harmonics  that  are  multiples  of  the  number  of  blades  are 
transmitted  (i.e.,  for  the  four-bladed  rotor,  only  0, 4,  8, 12,  etc.,  per  revolution  will  be  transmitted). 

Another  effect  to  be  noted  is  that  all  the  frequencies  are  coming  through  to  the  fixed  system  when 
a  blade  is  damaged  (i.e.,  1, 2, 3, 4, 5, 6/rev);  thus,  the  rotor  thrust  vector  is  tilting  at  these  frequencies 
and  shaking  the  helicopter.  In  Table  3,  a  basic  fundamental  regular  pattern  shows  up,  frequencies 
are  transmitted  into  the  fixed  system  plus  and  minus  1/rev  of  the  rotating  blade  forced  response 
frequency.  Also  the  blade  summing  process  yields  a  mechanical  frequency  filtering  effect  when  the 
blades  are  undamaged,  which  is  lost  when  a  blade  becomes  damaged.  An  explanation  of  the  BXak 
and  the  BYbk  terms  in  Tables  2  and  3,  where  k  =  1 . .  .6,  is  given  in  the  Appendix. 

5.  Pictorial  Representation  of  Rotor  Disk  Tilt 


Graphical  representations  of  the  fixed-system  BXa2  rotor  tilt  as  a  function  of  time  are  shown  in 
Figures  3-6  looking  in  the  negative  X  direction  into  the  Y-Z  plane. 
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ROTOR  DISK  PLANE 


Y  Projection 

Figure  3.  BXa2  Rotor  Disk  Plane  Tilt  at  t  =  0. 
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ROTOR  DISK  PLANE 


Y  Projection 

Figure  5.  BXa2  Rotor  Disk  Plane  Tilt  at  t  =  pi/(2  om). 
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Figure  6.  BXa2  Rotor  Disk  Plane  Tilt  at  t  =  3  pi/(4  om). 
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Figure  3  represents,  in  a  perspective  view,  a  picture  of  the  rotor  disk  that  a  viewer  would  see 
looking  straight  on  at  the  disk  from  the  X-axis  projected  on  the  Y-Z  plane.  This  shows  the 
undisturbed  disk  lying  horizontally,  and  also  the  tilted  disk  caused  by  a  BSa2  contribution  from  a 
30%  damaged  blade  response  to  1/rev  forcing.  This  picture  illustrates  the  magnitude  of  the  2/rev  tilt 
of  the  rotor  disk  in  the  fixed  system  at  time  =  0  due  to  a  damaged  blade’s  response  to  1/rev  forcing. 

In  Figure  4,  at  time  =  pi/(4  om),  the  disk  has  moved  back  down  to  the  undeflected  roll  tilt 
position.  Then  at  time  =  pi/(2  om),  the  rotor  disk  is  tilted  in  the  opposite  maximum  tilted  position. 
Finally,  at  time  =  3  pi/(4  om),  it  is  back  to  an  untilted  position  again.  Collectively,  Figures  3-6  show 
that  the  rotor  disk  will  go  through  a  complete  2/rev  cycle  in  pi/om  seconds,  or  half  a  rotor  cycle. 

From  Tables  2  and  3,  we  see  the  peculiarity  of  only  BXak  and  BYbk  coefficients  present,  with 
no  BXbk  nor  BYak  terms.  This  is  caused  by  blade  1  being  damaged;  if  blade  2  had  been  assumed 
damaged,  the  roles  would  have  been  reversed  (i.e.,  BXbk  and  BYak  terms  only  and  no  BXak  nor 
BYbk  terms).  See  the  Appendix  for  the  mathematical  details. 

6.  Summary 


An  elementary  helicopter  rotor  blade  dynamical  system  is  derived  to  demonstrate  the  elemental 
principles  by  which  blade  damage  affects  rotor  system  performance  and  can  affect  other  aspects, 
such  as  human  tolerance  to  vibrations,  structural  fatigue  life,  and  general  controllability. 

Starting  with  the  rotor’s  individual  blades  and  their  responses  to  forcing  (in  this  case,  harmonic 
sinusoidal  forces),  it  is  shown  they  have  an  integrated  effect  on  the  fixed-system  rotor  disk  tilt. 
When  the  blades  are  undamaged,  the  helicopter’s  trim  condition  is  determined  by  the  blade  1/rev 
flapping  of  similarly  responding  blades.  Also,  for  a  four-bladed  rotor,  there  is  an  inherent  natural 
4/rev  tilting  of  the  rotor  disk  in  the  fixed  system  that  is  always  present. 
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When  a  blade  is  damaged,  rotor  imbalance  exists  and  all  the  per  rev  frequencies  are  “loosed”  to 
be  transmitted  into  the  fixed  system.  A  plus  and  minus  1/rev  frequency  conversion  is  demonstrated 
on  the  rotor  disk  as  well  as  a  blade  summing  and  filtering  out  of  frequencies  for  an  undamaged  rotor. 

One  kind  of  blade  damage  is  studied,  that  of  the  loss  of  the  outer  portion  of  one  blade  of  a  set  of 
four.  Further  parametric  damages  can  be  imposed  on  the  model  to  determine  transmitted  vibration 
effects  as  a  function  of  damage  type. 

Since  a  helicopter  as  a  whole  responds  to  the  entire  rotor  rather  than  to  a  single  blade,  analysis 
of  the  helicopter  response  due  to  knowledge  of  single  blade  forcing  does  not  give  a  complete 
understanding  of  the  damaged  blade  effect  It  is  therefore  proposed  that  further  work  in  this  area  be 
undertaken. 
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Appendix: 

Explanation  of  BXak  and  BYbk  Coefficients 
Appearing  in  Tables  2  and  3 
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The  type  of  blade  response  from  this  analysis  is  of  the  form 


4  (  fi-n  (  (i-1) 

X  Cik  sink  w  +  2 n- -  sin  \|f  +  2n — - — 

i=i  4  4 


4  (  (i-1)  }  f  (i-1) 

X  Cy.  sink  \|t  +  2jt- - -  cos  t|r  +  271^— — 

i=i  4)1  4 


which  are  fixed-system  rotor  disk  roll  and  pitch  tilts  due  to  sin  k  t| r  rotor  blade  forcing. 


From  trigonometric  identifies 


(  (i-1)  ^  (  (i-1) 

X  CL  sink  ty  +  2k- - -  sin  \|/  +  2k — - — 

4  14 


1  f  (  (i-1)  >  (  (i-1) 

— X  Cy.  \ cos (k - 1)  Y  +  2n— —  -cos  (k  -1)  \\r  +  2k  — — 

2  I  4  l  4 


—  X  Cfc  cos(k-l)\|tcos  (k-1)  2k  ^  — 
2  4 


-X  C.±  sin (k  - 1)  \|f  sin  (k-l)27t^— 
2  4 


— X  0^  cos  (k  + 1)  \|t  cos  (k  + 1)  27t  — — —  + 
2  4 


-X  Cy.  sin  (k  + 1)  \|t  sin  (k  +  l)2:t^— ^ 
2  4 
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From  this  we  identify 


BXa(k-l)  =  —X  C^cosCk-l)  2k— — — 

2  4 

BXb(k-l)  =  --I  Cj,  sin (k - 1)  2k^^- 
2  4 

BXa(k  +  l)  =  X  C^cosCk  +  l)  2k— — — 

2  4 

BXb(k  +  l)  =  —  X  sin(k  +  l)  2n— — — 

2  4 

These  are  the  fixed-system  rotor  disk  rolling  terms.  In  a  similar  maimer,  the  pitching  terms  are 
calculated  as 


BYb (k  + 1)  =  -X  cos (k  + 1)  2n^-^ 

2  4 


BYa(k  +  l)  =  -X  Cj,  sin  (k  + 1)  2tt^— ^ 
2  4 


BYb(k-l)  =  -X  cos (k - 1)  2*«LJ> 

2  4 


BYa(k-l)  =  —X  C^sinCk-l)  2k— — — 
2  4 


Now,  for  instance,  for  k  =  2  (2/rev)  with  blade  1  damaged,  let  blade  1  have  an  amplitude  of  1.1, 
and  the  remaining  undamaged  blades  have  an  amplitude  of  1 .0.  Then  we  will  have 
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BXal  =  - 
2 


and 


1.1  cos(0)  +  1.0  cos 


it 


—  +1.0  cos  (it)  +  1.0  cos 
2 


=  i  {l.l  +  O-l.O  +0 }  =  |(0.1) 


BXbl  =  ~  jl.l  sin(0)  +  1.0  sin 


v2 ) 


+  1.0  sin  (  Itc)  +  1.0  sin 


f3n)' 


=  -i  {0  +  1.0  +  0  -  1.0 }=  0. 


For  blade  2  damage,  we  have 


BXal  =  -  {l.O  +  0-1.0  +  0}=0 


BXbl  =  --  {0  +  1.1  +  0  -  1.0}=  “(0.1). 
2  2 


We  see  here  how  the  damaged  blade  number  determines  whether  a  cosine  ak  or  sine  bk 
coefficient  is  obtained. 
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List  of  Terms 


Vi 

kz 

B  =  pi 

om  =  £2 
z 

Ml=Mi 
M2  =  M2 
Kz  =  Kz 
Kb  =  Kp 
e 
r 
F 
A; 

N 

t 

BX  =  Bx 
BY  =  By 
ak  =  ak 
bk  =  bk 
BXak 
BYbk 

Pi 

Q 


rotor  blade  azimuthal  displacement  of  ith  blade 

rotor  blade  harmonic  number 

ith  rotor  blade  flap  displacement 

rotor  angular  speed 

vertical  deflection  of  the  hub  mass 

rotor  blade  mass 

hub  mass 

hub  spring 

rotor  blade  spring 

rotor  blade  flap  hinge  offset 

rotor  blade  radius 

rotor  blade  forcing  function 

rotor  blade  forcing  amplitude  of  ith  blade 

number  of  blades  in  rotor 

time 

rotor  disk  lateral  tilt  in  the  fixed  system 
rotor  disk  longitudinal  tilt  in  the  fixed  system 
rotor  disk  cosine  coefficient  in  the  fixed  system 
rotor  disk  sine  coefficient  in  the  fixed  system 
rotor  disk  lateral  tilt  cosine  coefficient  kth  harmonic 
rotor  disk  longitudinal  tilt  sine  coefficient  kth  harmonic 

K 

amplitude  of  the  ith  rotor  blade  flap  displacement 
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